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ABSTRACT 

The  scattering  of  visible  light  in  hazes  and  fogs  was  investigated  experimen¬ 
tally,  using  balloon-borne  flash  bulb  clusters  as  sources  and  calibrated  cameras  as 
detectors.  By  defocusing  the  cameras,  measurements  of  both  the  direct  and  the 
scattered  components  could  be  obtained  on  the  same  film.  Several  detectors  were  sta¬ 
tioned  on  the  ground  at  various  distances  up  to  15  kms  from  the  source. 

Results  indicate  that  transmission  of  the  direct  components  decrease  exponen¬ 
tially  and  as  the  inverse  square  of  the  slant  range  R,  while  scattered  radiation  ti’ans- 
mission  can  be  described  by  a  different  exponential  and  the  inverse  first  power  of  the 
slant  range.  All  transmission  constants  can  be  related  to  the  visibility  at  the  time  of 

measurement.  Brightness  of  the  aureole  surrounding  the  source  decrease  inversely  as 
2  2 

R  sin  0  ,  where  ©  is  the  angular  distance  from  the  source  as  seen  by  the  observer. 
This  is  an  average  relationship  applying  at  all  distances  and  fog  densities  observed. 

A  general  expression  for  the  total  flux  from  a  point  source  reaching  an  observer 
is  given  in  terms  only  of  the  visibility  and  the  slant  range. 


I.  INTRODUCTION 


Evaluations  of  the  effect  of  various  meteorological  scattering  media  on  radiant 
energy  transport  in  the  atmosphere  have  been  attempted  both  theoretically  and  experi¬ 
mentally. 

Theoretical  studies  generally  follow  the  line  of  reasoning  which  specifies  model 
atmospheres  with  homogeneous  strata.  These  approaches  are  not  without  merit,  be¬ 
cause  they  indicate  the  manner  in  which  light  is  transmitted  through  various  scattering 
media.  The  weakness  in  these  approaches  seems  to  be  that  the  models  may  be  relatively 
unrealistic:  i.e.  ,  real  atmospheres  are  not  homogeneous,  and  radiant  energy  transfer 
may  be  influenced  by  losses  upward  and  out  of  the  atmosphere  and  by  reflections  from 
surfaces  whose  albedoes  are  not  well  known. 

Experimental  approaches  to  this  problem  can  be  instructive  by  indicating  the  flux 
density  received  from  a  point  source  through  distances  of  several  kilometers  for  those 
weather  situations  encountered.  The  diffuse  transmission  through  real  atmospheres  as 
measured  by  photometers  with  fields  of  view  of  about  150°  was  reported  earlier  by  the 
authors  (References  1  and  2).  This  study  indicated  the  general  trends  that  the  transmission 
assumes.  However,  the  relationship  between  the  direct  and  scattered  components  of  the 
total  transmission  is  unknown  because  the  photometers  measured  total  flux  only. 

The  subject  of  this  final  report  is  a  study  which  was  aimed  at  determining  the  be¬ 
havior  of  the  direct  and  scattered  components  of  the  total  transmission  through  various 
scattering  media,  as  well  as  the  possible  contribution  of  reflections  from  the  interven¬ 
ing  terrain.  The  results  are  discussed  in  terms  of  the  actual  measurements,  with  an  indi¬ 
cation  of  how  they  may  be  applied  to  direct  beam  measurements  made  in  similar  atmos¬ 
pheres  by  other  workers.  Finally  an  attempt  has  teen  made  to  relate  the  experimental 
results  to  theoretical  parameters  and  a  simple  model  atmosphere. 


-  1  - 


II.  INSTRUMENTATION 


The  principle  adopted  here  for  the  measurement  of  both  direct  and  scattered  com¬ 
ponents  of  the  flux  arriving  from  a  concentrated  source  was  one  in  which  both  could  be 
simultaneously  determined.  The  detection  equipment  consisted  simply  of  cameras  and 
film,  with  auxiliary  apparatus  needed  to  make  the  system  acceptably  quantitative.  Because 
the  light  sources  consisted  of  50  cm  diameter  clusters  of  photoflash  bulbs,  operated  at 
distances  of  at  least  several  hundred  meters,  their  images  would  have  been  extremely 
small  on  the  film  in  a  conventional  camera.  Moreover,  their  brightnesses  would  have  far 
exceeded  the  maximum  which  the  film  could  record.  To  obviate  these  difficulties,  the 
cameras  were  deliberately  defocused  to  the  point  where  the  source  produced  a  measurably 
large  circle  of  confusion  on  the  film.  The  effect  of  defocusing  on  the  resolution  of  the 
system  was  such  that  scattering  of  2°  or  less  could  not  be  distinguished  from  the  direct 
beam.  Since  the  wide  angle  scattering  was  of  major  interest,  this  loss  of  resolution  was 
judged  acceptable. 

Six  pairs  of  cameras  were  located  at  various  distances  from  a  balloon-borne  light 
source.  The  balloon  flight  was  tracked  with  the  camera  shutters  open  until  after  the  flash 
occurred.  As  a  result,  the  exposure  was  given  by  the  lens  aperture  times  the  flash  pulse 
duration,  superimposed  on  the  illuminance  of  the  background.  The  experiments  were  per¬ 
formed  at  night  to  reduce  the  background  illumination  to  a  minimum.  The  photographic 
technique  also  allowed  evaluation  of  the  background  brightness  and  other  terrestrial  lights. 

The  source  consisted  of  30  No.  50  photoflash  bulbs  symmetrically  arranged  and 
connected  to  fire  simultaneously.  Based  on  information  furnished  by  the  manufacturer, 
and  a  measured  configuration  efficiency  of  0.  85,  the  energy  emitted  by  the  source  was 
2.  4  x  10G  lumen-seconds  in  the  visible  spectrum,  or  12  kilocalories  over  the  entire  3800°K 
black-bodv  spectrum.  This  emerged  in  a  triangular  pulse  of  40  milliseconds  duration,  with 

g 

a  peak  emittance  of  1.  3  x  10  lumens,  or  2.4  megawatts. 

Films  exposed  during  the  field  experiments  were  previously  masked  along  the  edge. 
The  resulting  unexposed  section  of  film  was  later  exposed  to  a  calibrated  step  wedge  illu¬ 
minated  by  an  electronic  flash  source  whose  duration  had  been  made  equal  to  the  field  photo¬ 
flash  source.  These  films  were  then  developed  at  constant  temperature  with  nitrogen  burst 
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agitation.  All  negatives  exposed  to  the  same  source  were  developed  together  to  obtain  as 
much  uniformity  as  possible.  Analysis  of  each  sheet  of  film  consisted  of  reading  the  opti¬ 
cal  densities  of  the  step  wedge  and  constructing  an  H  and  D  curve.  The  measured  optical 
densities  of  the  flash  image  were  referred  to  this  H  and  D  curve  to  determine  the  corres¬ 
ponding  exposures. 

The  cameras  used  were  B  and  J  4  x  5  Press  View  Cameras  with  Schneider- 
Kreuznach  Xenar  f/3.5,  105  mm  lens  in  Synchro-Compur-P  shutter,  operated  about  1  cm 
off  focus.  The  film  pack  adapters  were  modified  to  mask  off  one  end  of  the  film  for  later 
exposure  to  the  calibrated  step  wedge.  All  exposures  in  the  field  were  made  at  f/3.  5,  with 
the  exposure  time  determined  by  the  pulse  duration  of  the  photoflash  bulbs. 

The  detector-source  combination  used,  Kodak  Tri-X  film  in  the  pack  form  and 
No.  50  photoflash  bulbs,  exhibit  relatively  flat  spectral  response  between  0.  42  and  0.  64 
micron,  with  a  slight  maximum  occurring  at  about  0.61  micron  wavelength.  The  film  re¬ 
sponse  cuts  off  at  0.  36  and  0.  67  micron  wavelength  (Reference  3). 

Beside  the  necessity  for  film  calibration,  use  of  cameras  and  film  for  photometry 
also  requires  determining  (1)  the  relationship  between  the  optical  density  and  the  source 
intensity;  (2)  the  efficiency  of  the  lens  as  a  function  of  incident  angle,  and  (3)  the  relative 
sensitivity  of  each  camera.  Each  of  these  factors  will  influence  the  accuracy  of  measure¬ 
ment  of  the  atmospheric  transmissions.  The  techniques  and  experiments  performed  to 
evaluate  and  control  these  problems  are  discussed  below. 

A.  FILM  CALIBRATION 

Calibration  of  Kodak  Tri-X  film  was  accomplished  by  exposing  an  edge  of  the  film 
to  a  gray  step  wedge  illuminated  by  a  modified  electronic  flash  unit.  The  electronic  flash 
tube  was  mounted  at  one  end  of  a  light-tight  blackened  tube  and  a  film  holder  was  attached 
at  the  other  end.  The  physical  features  of  the  flash  sensitometer  are  shown  in  Figure  1. 

A  sheet  of  film  is  placed  in  the  film  slot  so  that  the  masked  end  is  laid  against  the  calibra¬ 
ted  step  wedge.  The  felt-lined  cover  is  pressed  down,  holding  the  film  in  place,  and  the 
flash  tube  is  triggered.  The  film  is  then  removed  and  placed  on  a  hanger  for  eventual  de¬ 
velopment  along  with  other  films  in  that  particular  experimental  transmission  measure¬ 
ment  The  film  is  processed  in  total  darkness,  from  removal  from  the  film  pack  holder 
to  complete  "fixing"  after  development. 
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FELT  LINED  COVER 


FIGURE  I 

FLASH  SENSITOMETER 


Because  of  the  marked  difference  between  the  pulse  length  of  a  standard  electronic 
flash  and  that  of  a  photoflash  bulb ,  the  electronic  flash  was  modified  to  approximate  the 
time  duration  of  the  photoflash  bulbs.  By  adding  resistance  in  series  with  the  flash  lamp, 
the  discharge  was  slowed  down  to  about  30  milliseconds.  Independent  measurements  of 
the  pulse  length  in  the  visible  range  had  previously  indicated  this  to  be  the  approximate 
pulse  length  of  the  photoflash  bulbs. 

The  flash  sensitometer  was  tested  for  reproducibility  by  measuring  the  light  out¬ 
put  with  an  integrating  photometer.  The  variation  about  a  mean  value  was  2%  for  30 
measurements  of  the  total  visible  integrated  flux.  This  error  is  well  within  the  limits  to 
which  the  film  optical  densities  can  be  measured  with  the  densitometer  being  used. 

The  gray  step  wedge  used  in  the  sensitometer  was  a  standard  semi-transparent 
film  wedge  purchased  from  the  Welch  Scientific  Company.  This  wedge  is  normally  used 
to  calibrate  the  Welch  Densichron  instrument.  The  range  of  densities  was  increased  by 
covering  one  edge  of  the  wedge  with  a  neutral  density  filter  made  from  uniformly  exposed 
film.  The  features  of  the  modified  step  wedge  are  shown  in  Figure  2.  The  neutral  density 
filter  is  mounted  along  the  left  side  of  the  wedge;  the  numbers  on  the  right  edge  are  the 
step  wedge  numbers.  The  optical  densities,  P,  are  listed,  as  well  as  relative  log  exposure 
beside  each  step. 

All  negatives  from  which  data  were  derived  were  developed  with  nitrogen  burst 
agitation  at  a  constant  temperature.  An  indication  of  the  constancy  of  development  was 
obtained  by  monitoring  the  resultant  H  and  D  curves  for  each  sheet  of  film  as  a  function  of 
chronological  time.  The  films  were  developed  19  minutes  at  68°F  in  Kodak  DK  50  develope 
giving  a  gamma  of  about  0.  9.  A  typical  H  and  D  curve  is  shown  in  Figure  3.  The  relative 
log  exposure  of  each  step  (Figure  2)  is  indicated  on  the  abscissa  of  Figure  3.  The  overlap 
of  steps  is  a  consequence  of  the  choice  of  the  neutral  filter. 

A  comparison  oi  .12  H  and  D  curves  indicates  an  average  variation  of  0.026  optical 
density  unit  with  a  maximum  variation  of  0.  035  optical  density  unit  where  the  overlap 
occurs.  This  results  in  a  maximum  uncertainty  of  0.08  in  relative  log  exposure  units. 

Part  of  this  uncertainty  results  from  inaccuracy  in  measuring  the  optical  density,  and  pro¬ 
duces  a  variation  in  the  relative  log  exposure  of  0.05  unit.  This  value  is  considered  more 
representative,  because  the  overlap  magnifies  the  uncertainty.  The  atmospheric  trans- 
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NEUTRAL  STEP 
TENSITY  WEDGE 
FILTER  NUMBER 


Log  E 

4.89 

4.77 

4.71 

4.64 

4.51 
4.38 
4.24 
4.1  I 
3.95 

3.80 
3.61 
3.42 
3.23 
3.05 
2.88 
2.69 

2.52 
2.32 
2.16 
2.03 
1.90 


P 

0.07 
0.  16 
0.23 
0.32 
0.41 
0.57 
0.67 
0.80 
0.96 
1.12 
1.28 
1.45 
1.64 

1.81 
2.00 
2.19 
2.35 

2.53 
2.68 

2.81 
2.94 


FIGURE  2 

CALIBRATING  STEP  WEDGE 


mission  data  were  obtained  through  most  of  the  linear  portion  of  the  curve;  however,  be¬ 
cause  each  relative  log  exposure  is  derived  through  the  H  and  D  curve  for  each  sheet  of 
film,  data  can  be  obtained  throughout  the  H  and  D  curve  provided  the  optical  density  is 
above  the  fog  level  and  below  the  level  of  overexposure.  In  practice,  care  was  exercised 
to  determine  that  neither  of  these  limits  was  exceeded.  These  uncertainties  also  contri¬ 
bute  to  the  experimental  error,  and  are  reflected  in  the  atmosphere  transmissions. 

An  examination  was  made  of  the  effect  of  differences  in  color  temperatures  of  the 
electronic  flash  lamp  source  used  in  the  sensitometer  and  of  the  photoflash  bulb  sources 
used  in  the  field  experiments.  This  possibility  of  error  was  investigated  by  exposing  film 
to  both  sources  through  the  step  wedge.  Results  normalized  in  the  linear  portion  of  the  H 
and  D  curve  are  shown  diagramatically  in  Figure  4.  Except  at  the  lower  end  of  the  expo¬ 
sure  scale,  the  comparison  of  light  sources  indicates  that  the  use  of  an  electronic  flash 
lamp  source  in  the  sensitometer  introduces  only  a  negligible  error  in  the  calibrations  of 
the  film.  This  conclusion  is,  of  course,  valid  only  with  Kodak  Tri-X  (panchromatic)  film 
and  the  particular  electronic  flash  tube  and  No.  50  photoflash  bulbs  (color  temperatures 
about  6000°K  and  3800°K,  respectively). 

B.  CAMERA  CALIBRATION 

Calibration  of  the  cameras  required  the  measurement  of  two  major  parameters: 

(1)  lens  transmission  as  a  function  of  incident  angle  for  each  camera;  (2)  relative  sensiti¬ 
vity  of  each  camera. 

The  camera  lens  was  calibrated  by  making  photographs  of  a  fluorescent  lamp,  chosen 
because  of  its  relatively  uniform  brightness.  The  films  were  developed  with  nitrogen  burst 
agitation  after  exposure  to  the  calibration  wedge.  The  brightness  of  the  fluorescent  lamp 
was  measured  along  the  longitudinal  axis  at  the  points  whei'e  densitometric  measurements 
of  the  lamp's  image  were  to  be  made.  The  results  are  summarized  in  graphic  form  for  all 
cameras  in  Figure  5.  The  ordinate  is  the  lens  efficiency,  equal  to  the  normalized  ratio  of 
image  to  lamp  brightness;  the  abscissa  is  the  lateral  distance  from  the  center  of  the  film. 

It  may  be  seen  that  the  experimental  data  depart  from  the  theoretical  rate  of  decrease  pro¬ 
portional  to  cos^4>.  where  tan  4>  distance  from  center/image  distance  from  lens.  How¬ 
ever.  it  is  not  uncommon  for  experimental  lens  efficiencies  to  deviate  from  this  value 
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FIGURE  5 


(Reference  4).  The  experimental  curve  indicates  that  some  vignetting  is  present.  The 
experimental  lens  efficiency  curve  is  used  to  correct  the  source  image  density  accord¬ 
ing  to  the  image  distance  from  the  center  of  the  film. 

The  relative  sensitivity  of  each  camera  was  determined  by  observing  the  same 
flash  source  simultaneously  with  all  the  cameras  at  the  same  range  under  field  conditions. 
All  cameras  were  set  up  at  one  location  to  observe  seven  flashes  in  succession.  The  re¬ 
sulting  negatives  were  developed  as  described  and  the  relative  flux  density  for  each  flash 
was  obtained  by  the  procedure  detailed  below.  The  relative  sensitivity,  S,  for  any  one 
camera  is  defined  as  the  quotient  of  the  average  flux  density  measured  by  all  the  cameras 
and  the  flux  density  measured  by  the  given  camera.  For  each  camera,  standard  deviation 
of  the  individual  sensitivities  derived  from  each  flash  is  shown  in  Table  I  in  the  column 
headed  by  a.  .  This  is  a  measure  of  the  system  errors,  excluding  variations  in  the 
source  intensity,  and  it  can  be  seen  that  they  are  gratifyingly  small. 

TABLE  I 

CAMERA  RELATIVE  SENSITIVITIES 


Camera 

S 

ai  ’ 

1 

0.  96 

1.  23 

2 

1.05 

0.  80 

3 

1.  04 

0.60 

4 

1.  02 

0.57 

5 

1.01 

0.  57 

6 

0.95 

0.  93 

7 

0.99 

1. 06 

8 

0.96 

0.  90 

9 

0.  97 

0.  52 

10 

1.  02 

0.  97 

11 

1.01 

0.97 

12 

1.02 

1.  30 

Average 

1 . 00 

0.  87 
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C.  EXPERIMENTAL  PROCEDURE 


All  field  experiments  were  performed  at  night  on  Plum  Island,  Newburyport, 
Massachusetts.  (For  a  detailed  description  of  balloon  and  source  techniques  used 
here,  see  Reference  1. )  Cameras  were  mounted  in  pairs  on  a  swivel  support  so  that 
they  could  be  aimed  in  azimuth  and  elevation.  The  12  cameras  were  distributed  at 
five  or  six  different  distances  from  the  launching  site  of  the  balloon-borne  photoflash 
bulb  source.  When  the  balloon  was  released,  a  pilot  light  was  tracked  with  each  set 
of  cameras.  The  shutters  were  opened  after  launching  prior  to  the  flash,  and  closed 
immediately  after  the  flash.  The  balloon  was  also  tracked  with  two  theodolites  situa¬ 
ted  1.  26  km  apart.  Each  theodolite  determined  the  azimuth  and  elevation  angle  of  the 
flash  with  respect  to  the  mutual  base  line,  and  from  these  data  the  heights  of  the 
flashes  were  calculated  (see  Table  II).  Similarly,  the  horizontal  position  of  each  flash 
was  calculated  to  within  ±  0.  02  km,  and  from  this  the  slant  ranges  from  the  flash  to 
each  camera  were  determined  by  triangulation. 


TABLE  II 

MEASURED  HEIGHTS  OF  SOURCES 


Flash 

Height  (km) 

Deviation 

Atmosphere 

A 

0.  268 

±  0.013 

Light  haze 

B 

0.  225 

±  0. 015 

Light  haze 

C 

0.468 

A  0.023 

Light  haze 

D 

0.190 

A  0. 030 

Moderate  haze 

E 

0.  220 

A  0.015 

Dense  haze  and  smoke 

F 

0.19 

A  0.05 

Dense  fog 

G 

0. 14 

A  0.  08 

Dense  fog 

H 

0. 13 

t  0.  08 

Dense  fog 

The  heights  given  for  the  flash  in  dense  fog  are  less  accurate,  due  to  the  extreme 
difficulty  in  tracking.  The  pilot  light  was  visible  for  only  a  short  time  before  it  dis¬ 
appeared  into  the  fog,  and  indicated  only  roughly  where  the  flash  occured  in  space.  In 
those  cases  the  theodolite  was  used  to  determine  the  point  where  the  pilot  light  disappeared 
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from  view,  after  which  the  course  of  the  balloon  was  extrapolated  to  the  time  of  firing. 

The  largest  variable  in  the  experiment  was  the  brightness  of  the  photoflash  bulb 
source.  The  percentage  variance  of  the  source  brightness  was  determined  from  the  same 
data  previously  used  to  determine  the  individual  camera  sensitivities.  Using  these  seven 
flashes,  the  variance  of  the  flux  density  is  25%.  However,  because  the  flux  from  each 
flash  was  measured  at  several  ranges,  the  absolute  value  of  the  source  was  not  needed  to 
determine  relative  transmission  between  stations. 

When  each  field  experiment  was  completed,  the  films  from  each  camera  were  ex¬ 
tracted  from  the  film  packs  and  grouped  according  to  the  flash.  The  calibration  wedge 
was  exposed  on  each  film,  and  then  all  12  sheets  of  film  were  developed  together.  As  a 
result,  the  only  uncontrolled  variables  remaining  in  the  experiment  were  the  photoflash 
source  intensity  and  the  atmosphere.  Source  variations  could  be  eliminated  by  extra¬ 
polating  the  unscattered  component  to  zero  distance;  thus  the  measured  quantities  should 
be  functions  only  of  the  atmosphere. 
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III.  DATA  REDUCTION 


After  a  series  of  negatives  had  been  developed,  the  optical  density  of  the  central 
image,  ,  and  the  optical  density  at  specific  distances  from  the  central  image,  Pg  , 
were  measured  with  a  Welch  Densiehron  Densitometer.  The  general  features  of  negatives 
for  the  three  major  atmosphers  encountered  are  shown  in  Figure  6.  The  high  optical  den¬ 
sity  at  the  center  of  the  nearly  concentric  rings  represents  the  defocused  source  image; 
the  rings  about  this  image  represent  lines  of  equal  optical  density.  In  general,  the 
scattered  or  diffuse  energy  is  circularly  symmetric  about  the  source  in  atmospheres  of 
haze,  and  somewhat  elliptically  symmetrical  in  an  atmosphere  of  fog.  In  practice,  the 
optical  density  is  measured  in  the  four  quadrants  from  the  source  image  along  the  hori¬ 
zontal  and  vertical  axis  at  0.  25 ,  0.50,  0.75,  1.0,  1.5,  2.0,  3.0,  4.0,  etc.,  mm  from 
the  image  when  the  atmosphere  is  hazy.  However,  because  the  diffuse  energy  is  not 
symmetrical  in  foggy  atmospheres,  the  optical  density  of  a  negative  representing  a  source 
in  fog  is  measured  on  a  10-mm  grid  centered  on  the  source  image.  These  optical  densities 
are  converted  into  log  exposure  by  means  of  the  H  and  D  curve  for  that  particular  piece  of 
film.  The  unscattered  component  of  the  light  incident  on  the  camera  is  the  product  of  the 
exposure,  ED,  and  the  area  of  the  central  image,  A^.  The  scattered  energy  is  the  pro¬ 
duct  of  the  exposure,  Eg,  and  the  area  of  the  grid  of  density  measurements  observed. 

The  direct  and  scattered  fluxes,  respectively,  are  determined  as  follows: 


and 


where: 


D  .  —  E TA  S_F ; 

obs  D  D  C 


Sobs  “  ESASSCF; 


(3.1) 


(3.2) 


E^  and  E„  are  the  exposures  for  the  direct  and  scattered 
JL)  o 

energies  for  the  pulse  duration  of  the  photoflash  bulbs; 


A  is  the  area  over  which  the  optical  densities  and  p0 

D  S  are 

measured; 

S  is  the  relative  sensitivity  facl or  for  each  particular 
camera;  and 

F  is  the  filter  factor  when  a  neutral  density  filter  is  used  over 
the  lens. 
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The  neutral  density  filters  were  used  to  prevent  overexposure  when  the  source 
was  measured  at  ranges  less  than  1. 5  km.  These  filters  were  procured  from  Eastman 
Kodak  Company  who  specified  the  optical  densities  to  be  within  1%  of  the  values  1.  0, 

2.  0,  3.0  and  4.  0  for  each  filter. 

During  the  field  experimental  phase  of  this  study  10  flashes  were  fired  for 
atmospheric  transmission  measurements.  Two  of  these  were  discarded:  one  was  ob¬ 
served  by  an  insufficient  number  of  cameras  to  yield  significant  results;  the  other  was 
obviously  in  error  in  the  determination  of  the  flash  position  in  space.  Therefore,  eight 
measurements  of  the  total  transmission  are  presented:  three  for  very  light  haze,  one 
for  moderate  haze,  one  for  dense  haze  and  smoke,  and  three  for  dense  fog. 
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IV.  RANGE  DEPENDENCE  OF  DIRECT  AND  SCATTERED  FLUX 


Measurements  of  the  direct  and  scattered  flux  reaching  the  cameras  are  presented 

in  Figures  7  through  14  as  functions  of  the  slant  range,  R.  The  solid  curves  represent 

2 

the  data  multiplied  by  R  t.o  correct  for  inverse  square  dependence  and  were  thus  converted 
to  quantities  proportional  to  transmission.  The  direct  flux  in  all  cases  can  be  seen  to  be 
attenuated  exponentially  in  its  passage  through  the  atmosphere.  The  scattered  flux,  how¬ 
ever,  is  not  depleted  so  simply.  Theoretical  analysis  (Reference  5)  of  the  behavior  of 
radiation  singly  scattered  by  a  homogeneous  medium  indicates  that  this  intensity  should 
decrease  exponentially  and  approximately  inversely  with  the  distance,  especially  in  the 
neighborhood  of  the  source.  Based  on  this  theoretical  consideration,  the  scattered  com¬ 
ponents  were  multiplied  by  R  and  replotted  as  dashed  curves  on  Figures  7  through  14.  It 
is  clear  that  this  quantity  is  well  represented  by  a  straight  line  on  the  semi-logarithmic 
plots . 

Accordingly,  the  total  flux  arriving  at  a  detector  whose  slant  range  is  R  kilo¬ 
meters  from  an  isotropic  source  could  be  described  analytically  as  the  sum  of  two  compo¬ 
nents,  the  direct  flux  and  the  scattered  flux: 

T  =  D  +  S,  (4.1) 


where  T  is  the  total  flux,  D  and  S  are  the  direct  and  scattered  fluxes,  respectively .  Numer¬ 
ically,  D,  S,  and  T  are  the  fluxes  from  a  source  which  would  give  unit  flux  at  unit  distance 
in  vacuum.  The  direct  flux  is  then: 


D 


D  , 
obs 

(DR2)q 


e 


oR 


R 


2 


(4.2) 


where  cr  is  the  direct  flux  extinction  coefficient.  Results  for  each  test  are  normalized  by 

2 

dividing  all  the  observed  fluxes  by  the  extrapolated  value  of  (DR“)  at  zero  distance.  This 
appears  to  be  a  justifiable  procedure,  not  only  theoretically  but  also  as  a  result  of  Middle¬ 
ton's  experiments  (Reference  6);  it  results  from  the  assumption  of  strict  exponential  (and 
inverse  square)  attenuation  of  the  unscattered  flux  from  the  source.  Middleton's  work  also 
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shows  that  a  is  identical  with  the  atmospheric  extinction  coefficient. 


The  scattered  flux  may  be  described  as: 


S  = 


s  , 

obs 

(DR2)0 


(4.3) 


where  kg  is  the  scattered  flux  extinction  coefficient.  Combining  equations  (4.  2)  and 
(4.  3)  equation  (4. 1)  can  be  written  as: 


+ 


(4.4) 


One  other  quantity  of  interest  must  be  defined  in  the  same  terms.  The  apparent 
brightness  of  the  aureole  in  the  photographs  is  readily  derivable  from  the  densitometric 
data.  This  quantity,  B  (R,  9),  is  expressed  as  flux  per  steradian  subtended  at  the  detec¬ 
tor,  unit  flux  being  defined  as  before.  As  an  approximation,  the  brightness  was  averaged 
over  rings  of  constant  angle  0  away  from  the  line  joining  source  and  observer,  and  is 
generally  a  function  of  0  and  of  slant  distance,  R. 

The  experimental  values  of  a,  kg,  and  b  are  shown  in  Table  III  for  each  of  the 
eight  flashes.  These  quantities  were  obtained  by  a  least  squares  analysis  of  the  experimen¬ 
tal  data  shown  in  Figures  7  through  14,  with  the  exception  of  the  value  of  u  for  flashes  F 
and  G.  In  these  two  cases  the  visibility  was  so  low  that  only  the  stations  closest  to  the 

source  were  able  to  measure  the  direct  flux,  hence  the  unambiguous  straight  line  could  not 

2 

be  drawn  for  the  quantity  DR  vs  R  in  Figures  12  and  13.  As  an  approximation,  it  is 
assumed  that  the  source  intensities  for  flashes  F  and  G  were  identical  with  falsh  H,  taken 
within  two  hours  of  flash  F  and  one  hour  of  flash  G;  and  from  this  the  value  of  k^  and  b,  cr 
could  be  inferred. 
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FIGURE  8 
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FIGURE  12 
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TABLE  III 


TRANSMISSION  PARAMETERS 


Flash 

Meteorological 

Description 

a  (km-1) 

kg  (km  1) 

b 

A 

Light  haze 

0.  065 

0.  34 

0.  040 

B 

Light  haze 

0.043 

0.  38 

0.  085 

C 

Light  haze 

0.  052 

0.  44 

0.  047 

D 

Moderate  haze 

0. 16 

0.  43 

0.  009 

E 

Dense  haze  and  smoke 

0.58 

0.  96 

0.  52 

F 

Dense  fog 

15. 1 

1.  90 

1.4 

G 

Dense  fog 

13.  2 

1.91 

2.  2 

H 

Dense  fog 

15.  9 

2.  05 

3.  9 

It  was  gratifying  to  find  that  a  for  the  three  flashes  agreed  with  each  other  within  ±  10%. 

It  is  of  interest  to  determine  the  dependence  of  the  scattered  flux  extinction 
coefficient  upon  the  direct  flux  extinction  coefficient.  The  simplest  form  of  this  depend¬ 
ence,  as  illustrated  in  Figure  15,  is: 


As  the  parameter  b  of  equation  (4.  3)  is  quite  sensitive  to  the  experimental  data,  the  de¬ 
pendence  of  b  upon  the  direct  flux  extinction  coefficient  was  determined  for  those  atmos¬ 
pheres  which  obviously  had  significant  scattered  energy.  For  this  reason  the  dependence 
of  b  on  a  was  investigated  only  for  the  three  cases  of  fog,  dense  haze,  and  smoke,  and 
moderate  haze;  results  are  shown  in  Figure  16.  It  may  be  noted  that  this  does  not  intro¬ 
duce  excessive  error  for  the  three  cases  of  light  haze.  The  equation  of  a  simple  dependence 
of  the  parameter  b  upon  the  direct  flux  extinction  coefficient  is: 


FIGURE  16 


b  =  0. 6  ct 


1/2 


(4.6) 


Substitution  of  equations  (4.  5)  and  (4.  6)  into  equation  (4.  4)  yields  an  expression  for 
the  total  flux  from  an  isotropic  source  in  terms  of  the  atmospheric  extinction  coefficient 
and  the  distance: 


T 


-  ctR 


0.6  c 1/2  e  ~c1/3r 


R 


R 


(4.7) 


Equation  (4.  7)  is  empirical,  and  as  such  has  limitations  in  its  use.  Considering  its 
derivation  -  i.e.  ,  the  approximate  relationships  between  the  measured  characteristic  of 
the  eight  atmospheres  -  it  seems  reasonable  to  assume  that  this  equation  is  useful  for  atmos¬ 
pheric  extinction  coefficients  from  0.02  <  a  <20  km  . 

The  usefulness  of  equation  (4.  7)  for  computing  the  scattered  flux  is  illustrated  in 
Figure  17.  The  quantity  SR  vs  the  slant  range,  R,  is  plotted  for  the  experimental  data 
(reproduced  from  Figures  7  to  14)  at  the  left.  At  the  right,  the  computed  values  of  SR  are 
plotted  for  comparison.  It  can  readily  be  seen  that  the  computed  scattered  flux  approximates 
the  general  trends  quite  satisfactorily;  the  slopes  of  the  curve,  with  the  exception  of  the 
three  light  hazy  atmospheres,  are  well  reproduced.  The  magnitude  of  the  observed  scattered 
flux  in  haze  appears  to  be  low.  A  reasonable  explanation  is  that  this  flux  is  small,  and  as 
such  approaches  the  magnitude  of  the  background.  The  evaluation  of  the  scattered  flux  as  it 
approaches  the  magnitude  of  background  requires  more  than  the  simple  densitometric  tech¬ 
niques  used  in  this  preliminary  study.  As  the  optical  density  of  the  scattering  atmospheres 
increases,  there  is  a  commensurate  increase  in  the  amount  of  scattered  Ilux,  and  therefore 
less  apparent  deviation  between  the  experimental  and  calculated  scattered  fluxes. 

It  is  useful  to  compare  the  total  transmissions  predicted  by  equation  (4.  7)  with  the 
corresponding  observed  quantities  for  the  various  atmospheres.  The  total  energy  trans¬ 
missions  were  obtained  from  the  experimental  data  by  adding  the  direct  and  scattered  fluxes, 
DR2  and  SR2,  appearing  in  Figures  7  through  14.  Values  at  each  range  were  averaged  and 
plotted  on  Figure  18.  Curves  were  also  drawn  on  Figure  18,  representing  the  predictions 
of  equation  7  for  each  of  the  transmissions.  Agreement  between  experimental  and  predicted 
values  is  good. 
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FIGURE  17 


V.  ANGULAR  DEPENDENCE  OF  SCATTERED  ENERGY 


The  aureole  is  a  characteristic  phenomenon  when  light  sources  are  viewed  through 
hazy  or  foggy  atmospheres.  The  aureole  is,  of  course,  a  manifestation  of  light  scattered 
by  atmospheric  aerosols.  The  photographic  technique  described  earlier  in  this  report 
yields  information  on  the  angular  distribution  of  brightness  about  the  spherically  symme¬ 
tric  photoflash  sources  observed  under  various  conditions  of  atmospheric  optical  densities. 

It  is  the  aim  of  this  section  to  examine  the  angular  distribution  of  brightness  for  variability, 
for  significant  parameters,  and  for  consistency  (when  integrated,  this  brightness  function 
must  yield  the  total  scattered  flux. ) 

Theoretical  studies  of  the  angular  dependence  of  scattered  flux  have  been  based  on 
the  relationship  between  the  scatterer  and  the  wavelength  of  the  incident  radiation,  i.  e.  , 
on  the  Mie  theory.  Examples  of  this  approach  have  been  reported  by  Chu  and  Churchill, 
Lowan,  Van  de  Hulst  (References  7,  8  and  9)  and  others.  Unfortunately,  it  is  quite  diffi¬ 
cult  to  apply  this  theory  to  the  experimental  data  because  droplet  sizes  are  unknown,  and 
because  of  the  wide  range  of  atmospheric  optical  densities  measured.  For  this  reason  the 
approach  is  empirical. 

An  examination  of  the  angular  data  indicates  that  the  brightness  of  the  aureole 
varies  as  much  as  four  orders  of  magnitude  above  background.  A  large  amount  of  scatter 
was  observed  when  all  of  the  data  obtained  was  considered,  although  each  individual  photo¬ 
graphic  film  showed  a  regular  decrease  of  brightness  as  the  angle  from  the  source  increased. 

It  is  also  evident  that  if  the  distance  from  source  to  observer  is  maintained  constant, 
an  increase  in  the  optical  density  of  the  atmosphere  decreases  the  brightness  of  the  aureole 
near  the  source,  but  increases  the  angular  size  of  the  aureole.  This  trend  is  clearly  illus¬ 
trated  as  one  makes  the  transition  from  moderate  smoke  and  haze  to  fog.  It  is  not  apparent 
from  these  measurements  that  this  trend  is  continued  as  the  transition  is  made  from  mod¬ 
erate  smoke  and  haze  to  light  haze.  Alternate  explanations  for  this  change  in  trend  can  be 
sought;  the  first  being,  of  course,  that  for  light  haze  the  aureole  near  the  source  is,  in 
fact,  less  bright  than  the  aureole  at  the  same  angle  when  observed  in  moderate  haze.  This 
point  of  view  is  supported  if  one  postulates  that  the  particle  size  distribution  remains  con¬ 
stant  and  that  only  the  concentration  of  scatterers  increases  as  the  turbidity  increases  from 
light  haze  to  moderate  haze  and  smoke. 
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A  second  explanation  for  the  observed  increase  in  brightness  near  the  source  as 
the  optical  density  of  the  atmosphere  increases  (although  the  mean  free  path  is  still 
long)  attributes  this  increase  to  the  manner  in  which  the  brightness  measurements  were 
made.  The  densitometer  available  to  this  project  had  an  aperture  that  was  too  large  to 
yield  approximately  incremental  fluxes  when  measuring  aureoles  of  small  angular  extent. 
In  the  analysis  of  light  hazes,  where  the  brightness  falls  off  very  rapidly  with  angular 
distance,  the  average  flux  over  a  finite  incremental  angle  may  be  a  poor  approximation  to 
the  brightness  at  points  near  the  source.  This  effect,  plus  the  deliberate  defocussing  of 
the  cameras,  produced  an  average  flux  that  is  less  than  the  maximum  of  flux  per  unit 
solid  angle  received  over  the  incremental  angle. 

The  uncertainty  attendant  on  the  brightness  measurements  made  in  light  haze  has 
resulted  in  our  concentrating  attention  quantitatively  on  the  trends  from  moderate  dense 
haze  and  smoke,  to  fogs.  This  range  of  turbidity  is  the  region  where  the  aureole  is  pro¬ 
minent  and  changes  brightness  and  size  most  markedly.  It  is,  therefore,  the  region  of 
greatest  interest  both  in  theory  and  in  application. 

In  analyzing  the  data  of  brightness  as  a  function  of  observer  angle,  the  following 
line  of  reasoning  was  used.  Starting  from  the  subjective  impression  that  the  aureole 
appeared  as  a  luminous  mantle  surrounding  the  source,  particularly  when  viewed  at  an 
appreciable  distance,  the  brightness  data  taken  from  the  photographs  were  plotted  as  a 
function  of  the  apparent  lateral  distance  from  the  source.  This  distance  is  given  by  R 
sin  0,  where  R  is  the  source-to-observer  distance,  and  0  is  the  angle,  subtended  at  the 
observer,  between  the  direction  of  the  source  and  the  direction  in  which  brightness  was 
measured  (see  Figure  19).  Since  the  aureole  is  presumably  spherically  symmetrical, 
observation  at  a  given  value  of  R  sin  0  always  represents  the  same  path  through  it;  the 
only  difference  between  data  points  is  that  the  path  length  increases  with  increasing  dis¬ 
tance  from  the  source. 

The  empirical  relationship  between  brightness,  B,  and  the  lateral  distance  from 

2 

source,  R  sin  0,  is  shown  in  Figure  20,  where  log  B  is  plotted  against  log  (R  sin  0)  . 
Each  curve  represents  the  brightness  distribution  under  given  meteorological  conditions 
and  at  a  given  distance  from  the  source.  Most  of  the  curves  tend  toward  finite  maxima  at 
low  values  of  R  sin  0,  and  increase  in  uncertainty  at  larger  values  of  radial  distance  from 
the  source,  where  the  brightness  merges  into  the  background. 
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APPARENT  RADIAL  DISTANCE  FROM  SOURCE,  METERS 


It  can  be  seen  that  in  spite  of  the  fact  that  the  individual  curves  are  not  straight 
and  do  not  necessarily  approximate  a  45°  slope,  the  aggregation  follows  such  a  rela¬ 
tionship.  Within  a  factor  of  3,  the  bulk  of  the  data  can  be  represented  by  a  45°  line  in¬ 
dependent  of  atmospheric  conditions,  or,  analytically 

B  =  SL2Q4 -  flux  per  steradian  (5.1) 

R2  sin2© 

where  R  is  expressed  in  kilometers. 

The  fact  that  this  relationship  covers  a  range  of  over  four  decades  of  brightness 
(normalized  only  to  a  standard  source)  and  visibilities  varying  from  0.  3  to  100  km  is  in 
itself  a  remarkable  finding.  With  it,  one  can  calculate  the  amount  and  angular  distribu¬ 
tion  of  scattered  energy  transmitted  through  any  scattering  atmosphere  from  an  arbitray 
source  geometry,  and  obtain  an  answer  good  to  a  factor  of  three.  (It  should  be  noted 
that  the  direct  beam  contribution  must  be  computed  separately  to  obtain  the  illumination 
of  the  observer). 

This  expression  for  the  brightness  distribution  can,  of  course,  be  integrated  to 
obtain  the  total  scattered  flux  at  any  point.  From  the  definition  of  brightness,  the  flux 
received  on  a  unit  area  of  a  flat  detector  is  given  by: 

•  2 

sm  6^ 

S  =  \  B  cos  0  d  SI  =  7r  \  B  d  (sin20)  (5.2) 

J  Jsin  0 

o 


when  the  detector  face  is  oriented  normal  to  the  line  joining  source  and  detector,  and 
assuming  the  aureole  to  be  circularly  symmetrical  about  the  source.  Here  0o  represents 
the  minimum  angle  of  observation,  in  this  case  the  angle  set  by  the  defocusing  of  the 
cameras,  and  the  maximum  angle.  Because  the  brightness  decreases  rapidly  with 
angle,  this  maximum  may  be  assumed  to  be  sin  0^  =  1.  Introducing  the  emperical  bright¬ 
ness  relationship  and  integrating, 


S 


0. 004  ~ 
R2 


In 


(sin20o) 


(5.3) 
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From  the  known  camera  constants  and  the  amount  of  defocusing,  this  works  out  to  be: 

S  =  0.1/R2  (5.4) 

This  result,  is  at  variance  with  equation  (4.  3),  obtained  by  graphical  integration  of  the 
scattered  light  appearing  in  each  of  a  series  of  photographs  observing  a  given  flash. 
However,  within  the  error  of  a  factor  of  three  given  above,  the  equations  are  in  agree¬ 
ment  over  the  applicable  range .  An  idea  of  the  approximations  involved  may  be  gained  by 

2 

inspection  of  Figures  7  through  14;  equation  (5.4)  predicts  that  SR  =  0. 1,  which  is  very 
roughly  descriptive  of  the  course  of  the  data. 

Integration  of  the  angular  flux  distribution  thus  reveals  that  the  relationship  given 
by  equation  (5. 1)  is  only  an  approximation,  though  a  very  useful  one  for  giving  the  qualita¬ 
tive  characteristics  of  the  aureole  distributions.  It  is  believed  that  further  observations 
will  be  required  to  produce  more  refined  relationships.  However,  those  reported  here 
should  be  sufficient  to  permit  rough  calculations  to  be  made  of  the  flux  which  would  pene¬ 
trate  openings  of  particular  orientations  and  solid  angles. 


VI.  CONCLUSIONS 


The  preliminary  atmospheric  transmissions  determined  with  the  photographic 
technique  were  made  in  atmospheres  covering  a  wide  range  of  visibilities.  From  these 
eight  transmissions  it  has  been  possible  to  derive  an  empirical  relationship  between  the 
scattered  and  direct  flux  extinction  coefficient.  From  this,  an  empirical  equation  for  the 
total  transmission  from  an  isotropic  source  may  be  computed,  as  it  was  shown  in  Section 
IV,  equation  (4.7): 


a  R 


T  = 


R 


0.6,t1/2  e  -lrl/3  R 
R 


(4.7) 


In  equation  (4.  7)  the  two  independent  variables  are  the  range  in  km  and  the  optical  density 
of  the  intervening  atmosphere.  The  optical  density  of  the  atmosphere  is  expressed  as  the 
extinction  coefficient.  Visibility  in  km  may  be  substituted  for  extinction  coefficient,  a  , 
using  the  standard  relationship  for  the  visual  range: 


V 

m 


3.912 

<7 


(6.1) 


The  standard  value  of  0.  02  is  used  for  C,  the  threshold  of  brightness  contrast. 

The  angular  distribution  of  the  scattered  flux  arriving  at  the  observer  was  likewise 
derived  from  the  photographs.  It  was  found  that  the  apparent  brightness  seen  by  the  observer 
could  be  approximately  related  to  the  apparent  radial  distance  R  (in  km)  from  the  source: 
within  a  factor  of  three,  this  brightness  could  be  given  by: 

g  -  flux  per  ste radian  (6-2) 

r 

independently  of  the  atmospheric  transmission. 

The  empirical  expressions  constructed  lor  the  direct,  total  scattered,  and  angular 
distribution  of  flux  lead  to  precise  statements  of  the  conditions  under  which  the  scattered 
flux  predominates  over  the  direct  flux.  For  atmospheres  where  the  visibility  is  less  than 
•1  km  (i.e.  a  >  1  km  *),  the  scattered  flux  becomes  equal  to  the  flux  in  the  direct  beam 
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somewhere  at  a  range  that  is  less  than  1.  5  km.  It  is  only  under  conditions  of  good  visibil¬ 
ity  that  scattered  flux  is  not  a  significant  fraction  of  the  total  radiant  energy.  As  a  rule  of 
thumb,  the  direct  beam  contributes  the  main  fraction  of  the  radiant  energy  when  visibilities 
are  greater  than  40  km  (  <r  =  0.  1  km  )  while  flux  contribution  from  the  direct  beam  may  be 
neglected  in  respect  to  the  contribution  from  the  scattered  flux  for  visibilities  less  than  0.  4 
km  (  a  =  10  km  ^).  For  intermediate  visibilities,  both  the  direct  and  scattered  flux  must  be 
considered  in  estimates  of  the  total  flux.  In  this  range  of  visibility,  the  fraction  of  each  flux 
that  will  be  received  by  a  given  detector  depends  not  only  on  the  extinction  properties  of  the 
atmosphere  (i.  e,  cr  dependent)  but  also  on  the  distance  from  the  source  at  which  the  measure¬ 
ment  is  made  (i.e.  ,  R  dependent). 

As  an  illustration  of  the  <7  and  R  dependence,  consider  what  type  of  atmosphere 
scatters  the  greatest  amount  of  energy  to  a  fixed  detector.  This  question  can  be  answered 
by  considering  equation  (4.  7).  Differentiating  the  second  term  with  respect  to  cr,  keeping  R 
constant,  and  setting  the  resulting  expression  equal  to  zero,  the  functional  relationship  be¬ 
comes: 

^  =  ^2R  <6'3> 

Equation  (6.  3)  states  that  the  atmosphere  scattering  the  greatest  amount  of  flux  to  a  detector 
depends  not  only  on  the  extinction  properties  of  the  atmosphere,  but  also  on  the  distance  be¬ 
tween  detector  and  source.  For  example,  a  detector  1  km  from  the  source  would  receive  the 
greatest  amount  of  scattered  flux  from  an  atmosphere  where  a  =  3.4  km  .  Range,  in  addi¬ 
tion  to  the  extinction  coefficient,  enters  into  a  proper  description  of  the  properties  of  a 
scattering  atmosphere. 

Enough  work  was  done  with  the  photographic  technique  to  establish  its  utility  as  a 
method  for  studying  light  scattering  in  the  atmosphere  and  to  show  where  it  needed  improve¬ 
ment  to  yield  more  precise  results.  It  is  believed  that  further  measurements  with  brighter 
sources,  more  accurate  distance  determinations,  and  better  densitometry  equipment  will 
permit  firmer  conclusions  to  be  drawn  about  the  scattering  properties  of  real  atmospheres. 
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